Quantum walks with memory — goldfish, elephants and wise old men 
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Quantum walks have emerged as an interesting approach to quantum information processing, 
exhibiting many unique properties compared to the analogous classical random walk. Here we in- 
troduce a model for a discrete-time quantum walk with memory by endowing the walker with 
multiple recycled coins and using a physical memory function via a history dependent coin flip. By 
numerical simulation we observe several phenomena. First in one dimension, walkers with memory 
have persistent quantum ballistic speed up over classical walks just as found in previous studies of 
multi-coined walks with trivial memory function. However, measurement of the multi-coin state can 
dramatically shift the mean of the spatial distribution. Second, we consider spatial entanglement in a 
two-dimensional quantum walk with memory and find that memory destroys entanglement between 
the spatial dimensions, even when entangling coins are employed. Finally, we explore behaviour in 
the presence of spatial randomness and find that in contrast to single coined walks, multi-coined 
walks do not localise and in fact a memory function can speed up the walk relative to a fully de- 
cohered multi-coin walker with trivial memory. We explicitly show how to construct linear optics 
circuits implementing the walks, and discuss prospects for classical simulation. 



I. INTRODUCTION 

Quantum computation [ ] is believed to allow cer- 
tain computational problems to be solved much more 
quickly than on classical computers. The standard ap- 
proach to quantum computation is using the circuit 
model, whereby an algorithm is decomposed into a set of 
quantum gates. Another contender is the measurement 
based model [ , ], whereby a highly entangled state is 
prepared and computation proceeds via only single qubit 
measurements. More recently, the quantum walk formal- 
ism, the quantum analogy of a classical random walk, 
has emerged as an interesting approach to implement- 
ing quantum computational tasks [4-7]. Here a walker^ 
i.e. a particle, is located at vertices in a graph and is al- 
lowed to coherently 'hop' along the edges of the graph 
to other vertices. This approach has proved fruitful for 
algorithm design [8-16] and is known to be universal for 
quantum computation [ , ]. Numerous experiments 
have begun to demonstrate elementary quantum walks, 
particularly photonic implementations where the walker 
is a single photon [19-26]. The standard models consider 
single walkers. However steps have been made towards 
implementing multi- walker quantum walks [22, 27, 28]. 

In the usual discrete-time quantum walk model a 
walker has a coin value associated with it, which speci- 
fies the direction the walker will take when propagating 
through the graph. In this paper we consider the case 
where multiple coins are employed, which can be inter- 
preted as memory of previous coin values. We show that 
this memory drastically affects the evolution of the quan- 
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tum walk, and indeed can result in a transition from 
quantum statistics to classical statistics. We study the 
effects of history-dependent coins and find that this can 
lead to various diffusive phenomena. We consider the ef- 
fect that measuring memory elements has on the evo- 
lution of the walk, and discuss the rate of spread of 
the probability distribution as a function of the memory 
length. We explicitly show how to experimentally con- 
struct an optical quantum walk with memory and discuss 
the challenges of classical simulation of such systems. 

Numerous authors have considered memory effects in 
the context of classical random walks [29-37], and some 
steps have also been made in the quantum context [38- 

], which we build upon. 



II. QUANTUM WALK FORMALISM FOR 
LINEAR GRAPHS 

We begin by defining a walker as a bipartite system, 
|x,c), where x is the position of the walker in a graph 
and c is the coin value which dictates the direction of 
the walker. In the standard discrete-time quantum walk 
formalism for a single walker on a linear graph, the evo- 
lution is decomposed into two steps, coin (C) and step 
(iS), defined as, 

C: |x,c)^5]45|x,i), 

3 

S : \x,c) ^ |x + c, c), (1) 

where A^^^ is a unitary coin matrix defining the transition 
amplitudes at position x. The coin value takes values ±1 
(right or left respectively). The coin operator coherently 
manipulates the coin value, putting it into a superposi- 
tion of left and right, while the step operator updates the 
position value according to the newly chosen coin value. 
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After the coin operator, the coin value can be interpreted 
as the direction the walker will propagate at the the next 
step, whereas after the step operator it can be regarded 
as memory of the previous coin direction. After evolv- 
ing t time steps the output state is |V^out) = (5'C)^|?/^in). 
This formalism could be logically extended to operate on 
graphs of higher order by extending the coin space be- 
yond ±1, but we will initially focus on the case of linear 
graphs for ease of exposition. 

The rate of spread of a quantum walk is quantified by 
the variance: 

aHt) = ^Mt){i-fir. (2) 

where /i = ^^Pii^ and Pi{t) is the probability to measure 
the walker at position i after a number of time steps t. 
A classical random walk on a line has variance cr^(t) = t 
whereas the quantum walk on a line cr^(t) oc where 
the proportionality constant depends on the initial state 
and the choice of coin rotation [ ]. In the presence of 
decoherence, either acting on the coin or spatial degrees 
of freedom, the quantum walk eventually behaves like a 
classical random walk with linear dispersion [ 1]. 

III. FORMALISM WITH MEMORY 

In the case of an arbitrary number of memory elements 
we introduce additional memory Hilbert spaces. We fo- 
cus on discrete-time quantum walks since no coins are 
present in the continuous-time quantum walk, and there- 
fore memory does not naturally arise. With N memory 
elements the walker is defined as |x, ci, . . . , cat), and we 
decompose the evolution into three stages by adding an 
additional memory update (M) operator. Then the evo- 
lution can be decomposed into, 

C : |x,Ci,...,CAr) ^A^^^^^.|x,Ci,...,CAr-i,j), 

3 

S \ |x, ci, . . . , Cat) ^ |x + Cat, Ci, . . . , cat), 
M : |x,ci, . . . ,CAr) ^ |x,CAr,Ci, . . . ,CAr-i), (3) 

where q is the coin value at memory element z, i.e. mem- 
ory of the coin value i steps ago. All the operations can 
be done quasi-locally in space (see Fig. 1). The coin op- 
erator coherently manipulates the last memory element 
into a superposition of left and right, the step opera- 
tor updates the current position according to the last 
memory element, and the memory update operator cycli- 
cally permutes the memory elements such that q q+i 
(i.e. what was the iih memory element next becomes the 
(i + l)th memory element). When = 1 we refer to this 
as a goldfish walk (goldfish have poor short-term mem- 
ory) , and when N = twe refer to this as an elephant walk 
(elephants have exceptional long-term memory). 

Our definition of a walk with memory is similar to that 
by Brun et al. [38] who introduced the use of cyclic per- 
mutations to update the memory registers but without a 




FIG. 1: Schematic of the quantum walk with memory. The 
walker can be viewed as a molecule carrying N spin- 1/2 par- 
ticles. The coin flip operation on spin A" is a unitary rotation 
with angle that is determined by the total polarisation of the 
other A — 1 memory spins (i.e. a memory function). As de- 
scribed in Sec. VI A the unitary can be generated by an Ising 
interaction with the indicated coupling graph. The centre spin 
N then performs the conditional shift operation on the entire 
molecule in space. Finally, in the memory update step, the 
spins are cyclically permuted. 

memory dependent coin flip. Flitney et al. [ ], consid- 
ered the case of some form of memory dependent coin flip 
in the context of Parrondo games. See also McGettrick et 
al. [40] for a treatment of quantum walks with a length- 
two memory. One of the key results of that work is that 
quantum walks with recycled coins do not decohere pro- 
vided that the coins are not recycled faster than every 
other time step. If, however, the number of coins does 
equal the number of time steps then the walks converges 
to the classical random walk behaviour. Here we give the 
first systematic numerical study of the effect of memory 
dependent coin flips on the behaviour of quantum walks. 

One might ask why we have chosen the coin operator 
to manipulate the last memory element. This is chosen 
to enforce the interpretation that q represents memory 
of the coin i steps ago. If, say, the first memory element 
were updated, this interpretation would no longer apply. 

After evolving t time steps the output state is 
I'^out) = {MSCy\ilJin)- To enforce elegant symmetry into 
the output state henceforth we symmetrize the input 
state to be of the form, 

IV-in} = |0) ® (I - 1, . . . , -1) + I + 1, . . . , +1»/V2, (4) 
and additionally employ a balanced coin, 

A(^) = e-'f = -L ( T) Va:, (5) 

which applies equal transition amplitudes in the left and 
right directions at each position. This particular form of 
highly entangled input for the coin state is convenient 
for visualisation purposes but is not necessary. In fact a 
product state input of the form |0) (8) | — 1)*^^ produces 
qualitatively similar behaviour with respect to the vari- 
ance. 

The full time evolution over 12 steps of the limiting 
goldfish and elephant walks are shown in Fig. 2. The final 
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output position probability distributions (summed over 
all coin values) for different memory sizes are illustrated 
in Fig. 3. In the case of a goldfish walk we observe the 
usual double-peaked distribution of a standard quantum 
walk with a single coin, and in the elephant walk we 
observe the classical statistics of a binomial distribution. 

The reason we observe classical statistics for the ele- 
phant walk is that the walker always has full memory 
of its entire evolution. Thus, no two trajectories inter- 
fere and therefore evolve independently of one another, 
giving rise to a classical distribution. An alternate way 
of thinking about this is that for t < N we are effec- 
tively using a fresh coin at each step. Thus, we observe 
classical behaviour for t < after which quantum be- 
haviour emerges. In general there is a high level of entan- 
glement between the position and memory subsystems. 
Note the elephant walk yields the same distribution one 
would obtain for a goldfish walk where the coin value 
is measured at each time step, or equivalent ly the coin 
value undergoes decoherence at each time step, as was 
demonstrated experimentally by Broome et al. [ ]. Im- 
portantly, even though we observe classical statistics in 
the elephant walk, our state is nonetheless a highly en- 
tangled pure state and no decoherence has taken place. 
The evolution is unitary and by reversing time we can al- 
ways evolve back to the input state, which is not possible 
in the truly classical case. 




IV. COIN MEASUREMENT & 
DECOHERENECE 

We now consider the evolution of the quantum walk 
when we measure some of the coins. In Fig. 4 we plot 
the output distribution for a an elephant walk with and 
without measurement of the last four coins in the mem- 
ory. After measurement we are left with a usual quantum 
walk distribution, which has been shifted according to 
the measurement outcomes of the memory registers and 
allowed to run for a shorter period of time. Specifically, 
measuring the last k memory registers, with outcomes 
Oi 'resets' the quantum walk and shifts its evolution to 
have an origin at = X^iLi allowed to run for 

t — k steps. This behaviour would ordinarily be observed 
in a standard quantum walk if measurements were per- 
formed intermittently during the evolution. But here the 
measurement can be performed after the evolution, which 
projects the distribution into a region constrained by the 
measurement outcomes. In the case of an elephant walk, 
if we measure all the coins we can effectively reconstruct 
the entire path the walker followed to reach its destina- 
tion. 

If decoherence is present in the memory registers, this 
is equivalent to the environment measuring the respective 
registers without giving us information about the mea- 
surement outcome. Thus, decoherence in the memory will 
yield a mixed state over the distributions obtained by dif- 
ferent memory measurement outcomes. This can lead to 
coherence over different length scales. Suppose the envi- 
ronment measures the last coins and is able to distinguish 
between (1, . . . , 1) and ( — 1, . . . , —1), which we acknowl- 
edge is a little artificial. This will yield a mixture of two 
coherent distributions located at different origins on the 
line, as per Fig. 4. Each sub-distribution exhibits full co- 
herence, whereas the total distribution is a highly mixed 
state. If the two sub-distributions are non-overlapping, 
tracing out the positions associated with one distribu- 
tion will yield a pure state over the remaining positions. 

See Refs. [42, 43] for a discussion on quantum walks 
with decoherent coins, and Broome et al. [ ] for an ex- 
perimental photonic implementation of a quantum walk 
with decoherent coins showing the transition from quan- 
tum to classical as the decoherence is increased. 



FIG. 2: (Colour online) Time evolution over 12 steps in the 
cases of a goldfish walk (N = 1), where we observe the usual 
double-peaked quantum walk ballistic spreading (top) and an 
elephant walk (N — t — 12), where we observe classical ran- 
dom walk statistics (bottom) . Note that even though classical 
statistics are observed in the elephant walk, no decoherence 
has taken place and the state is still pure. This is the limiting 
case of the Brun et al. result for walks with an equal number 
of coins and walk steps. 



V. TWO SPATIAL DIMENSIONS 

Next we turn our attention to a quantum walk on a two 
dimensional lattice. This is easily constructed by taking 
the formalism introduced earlier and expanding the coin 
degree of freedom to have four basis states (up, down, 
left and right). Now the basis states are of the form 
|x, ci, . . . , Cat), where x and y denote the two spatial 
dimensions. We begin with the walker localised at the 
centre of the graph and consider the cases where the coin 
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FIG. 3: (Colour online) Output distributions for t = 12, with different memory lengths N. From left to right: N = 1 (goldfish 
walk), = 3, = 10 and = t = 12 (elephant walk). In the case of a goldfish walk we observe the usual double-peaked 
distribution of a standard single-coined quantum walk, and in the case of an elephant walk we observe the classical statistics 
of a binomial distribution. 



FIG. 4: (Colour online) Output distributions for an elephant walk (N — t — 12). (left) no measurement, (middle) measurement 
of the last four coins with outcomes (1, 1, 1, 1), (right) measurement of the last four coins with outcomes ( — 1,-1,— 1,-1). 



is separable 



A 



(x) _ 
2D ~ 



(6) 



where e~*^^ is the coin from Eq. 5 (i.e. the coin acts on 
the two spatial dimensions independently) or maximally 
entangling, 



^2D ent ~ 



/ 1 \ 

10 

1 

\0 1 0/ 



Vx, (7) 



both with and without memory. The final probability dis- 
tributions are illustrated in Fig. 5. It can be visually seen 
that all the distributions are separable across the two 
spatial degrees of freedom, except in the case where en- 
tangling coins are employed and there is no memory. 

Interestingly, even when entangling coins are em- 
ployed, we do not observe spatial entanglement when 
memory is introduced. Thus, not only does memory re- 
duce the walk to a classical binomial distribution, but it 
also destroys entanglement between the two spatial de- 
grees of freedom. When the probability distribution is 
separable this implies the walker's x value is indepen- 
dent of its y value and we will always observe the same 
marginal x distribution. Thus, measurement along one di- 
rection has no impact on the measured distribution in the 
other. On the other hand, with spatial entanglement (i.e. 
inseparability), measurement in the x direction projects 
the y distribution onto a function of the measured x out- 
come. 

To quantify the spatial entanglement, in Fig. 6 we plot 
the spatial entanglement dynamics against time, con- 
firming that spatial entanglement only persists with en- 



tangling coins and no memory. The entanglement met- 
ric used is that described in Ref. [27]. We diagonalise 
the spatial probability distribution matrix and then cal- 
culate the Shannon entropy of the diagonal elements, 
S = — Xi log2 Xi. When the spatial distribution is sep- 
arable (i.e. there is no spatial entanglement), there is only 
one diagonal element and 5* = 0. For non-separable dis- 
tributions we have 5 > 0, a signature of entanglement. 



VI. HISTORY-DEPENDENT COINS 

Until now we have considered coins that are indepen- 
dent of the full memory history and only depend on 
the last memory element. We now consider the scenario 
where the coin operator is a function of the full previous 
history. We refer to this as a wise old man walk (wise old 
men have varying degrees of memory - some good, some 
horrendous - but they always base their actions on their 
past). 

Specifically, we will assume biased coins of the form. 



A(^)(i9) 



COS sin 
— sin^ cos^ 



(8) 



where the angle is computed from the memory history 

by 



TT 
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4 N-1 



(9) 



and < (j) < 1 is a parameter which specifies to what 
extent the memory is taken into account when choosing 
the coin operator. We refer to as the memory function, 
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FIG. 5: (Colour online) Spatial distribution of a quantum walk on a 2D lattice with t = 6 using separable (Eq. 6) and entangling 
(Eq. 7) coin flips: (1) Goldfish walk, separable coins; (2) Elephant walk, separable coins; (3) Goldfish walk, entangling coins; (4) 
Elephant walk, entangling coins. The distribution is separable across the two spatial degrees of freedom when either separable 
coins are employed or there is full memory. The only instance where the distribution is non-separable (i.e. where there is 
entanglement between the two spatial degrees of freedom) is when entangling coins are employed and there is no memory. 






Separable goldfish 




Entangling goldfish 




Separable elephant 




Entangling elephant 



FIG. 6: (Colour online) Spatial entanglement dynamics of 2D 
goldfish and elephant quantum walks, with separable and en- 
tangling coins. We only observe spatial entanglement when 
the coin is entangling and there is no memory. 



and it determines the diffusive properties of the walk. 
Of course this definition only works when N > 2. When 
^ = the coin operator reduces to a balanced Hadamard 
coin, independent of the history, and we reduce to walks 
similar to those previously discussed. 

In Fig. 7 we plot the time evolution of a wise old man 
walk with N = 5 for different history parameters ^ and 
memory function from Eq. 9. This time we use the in- 
put state |?/^in) = |0, 1, . . . , 1), i.e. we do not symmetrize 
the input state . When ^ = we reduce to a normal 
Hadamard walk. As (j) increases we observe reduced dis- 
persion, biased in one direction, yielding an asymmetric 
position distribution. To quantify this, in Fig. 8 we plot 
the variance of the distribution against time for differ- 
ent values of ^, confirming that higher <p results in re- 
duced dispersion. All distributions in Fig. 8 closely fit to 
quadratic curves. 

Fig. 7 can be understood as follows. Assume ^ = 1. 
If the walker has a long history of moving to the right, 
6 7r/2, and the coin operator reduces to the bit-flip 

matrix, A^^\^) = X = J J 

favours switching direction and heading back towards the 
origin. On the other hand, if the walker has a long history 
of moving to the left, ^ ^ 0, the coin operator reduces 
to the identity matrix, A^^\0) = I = diag(l, 1) Vx, and 
the walker will favour continuing its present course. 



Vx. Then the walker 



Thus in one direction the walker resists further dif- 
fusion, whilst in the other it does not, giving rise to 
the asymmetric distribution seen in Fig. 7. The latter 
property is better understood by considering the input 
state |?/^in) = |0, —1, . . . , —1). In this case the coin is the 
identity matrix, and remains so, and the walker lin- 
early shoots to the left without entering a superposition 
(graphic not shown). For this reason in the Fig. 7 we 
have only considered the \^|Jin) = term, which 

nicely illustrates the resistive diffusive behaviour in the 
rightward direction. 

Of course there is nothing unique about our choice of 
memory function in Eq. 9. Any function could be chosen, 
giving rise to a plethora of different diffusive phenomena. 
Another simple example is to choose the memory func- 
tion to be, 
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TV- 



(10) 



In this instance the coin reduces to a bit-flip only when 
the walker's history is consistently in the same direc- 
tion and to a Hadamard when the history is balanced. 
Thus, this memory function will resist excessive diffusion 
in both directions, yielding the much more localised dis- 
tribution shown in Fig. 9 (once again with a symmetrized 
input state). 

The two memory functions discussed previously had 
the effect of reducing dispersion. One can also choose 
memory functions which enhance diffusion. If we choose 
our memory function to be. 



TT 

4 



4 N-1 



(11) 



diffusion will be enhanced for larger cj). When the history 
is consistently in the same direction the coin reduces to 
the identity, and to a Hadamard when the history is bal- 
anced. This is shown in Fig. 10 (with a symmetrized input 
state). In Fig. 11 we plot the variance against time for 
various with memory function from Eq. 11. Once again 
all curves fit closely to quadratic distributions. 

The memory function from Eq. 11 is just Eq. 10 with 
a negative value of (j). Thus, for the memory function in 



6 



N-5,0=O^ N-5,(^^l 




Positiciri Posilion Posiliori 



FIG. 7: (Colour online) Time evolution over 12 steps for a wise old man walk with N = 5 and memory function from Eq. 9. 
(left) (/) = (a Hadamard walk), (middle) = 0.5 and (right) = 1. As increases, diffusion in one direction becomes more 
limited than the other, yielding an asymmetric distribution. 
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FIG. 8: (Colour online) Evolution of the variance of the prob- 
ability distribution against time for Suu N — 5 wise old man 
walk for different history parameters (j) and memory function 
from Eq. 9. As (f) increases, diffusion decreases. All curves fit 
closely to quadratic distributions. 
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FIG. 9: Time evolution over 12 steps for a wise old man walk 
with N — 5 and memory function from Eq. 10. Diffusion is 
restricted in both directions, yielding a much more localised 
distribution than Fig. 7(left). 



Eq. 10 there are three regimes. When ^ = we observe 
usual Hadamard ballistic spreading. When < < 1 we 
observe reduced diffusion. And when — 1 < ^ < we ob- 
serve enhanced diffusion. Note that cx V (j). 



FIG. 10: Time evolution over 12 steps for a wise old man 
walk with N = 5 and memory function from Eq. 11. (top) 
(f) — 0.5, (bottom) (f) — 1. Diffusion is enhanced in both di- 
rections compared to Fig. 7(left). For = 1 we observe that 
the walker linearly shoots off in both directions. In general 
this will not be the case, but comes as a result of the choice 
of initial state |^in) = (|0, -1, . . . , -1) + |0, +1, . . . , +1))/^, 
yielding the maximum possible dispersion. 



A. A physical model for history dependent coined 
walks 

Let us consider how to physically realise a history de- 
pendent coin rotation operation. For the memory func- 
tion defined in Eq. 9 we need to perform a conditional 
rotation on the Nth coin based on the states of the N — 1 
coins that specify the remembered past of length — 1 ^ . 



^ Physical realisations of memory functions Eqs. 10, 11 are more 
difficult to achieve and we do not discuss them here. 
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N=5 
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6=0.4 




6=0.6 




6=0.8 




6=1.0 



FIG. 11: (Colour online) Evolution of the variance of the prob- 
ability distribution against time for an N — 5 wise old man 
walk for different history parameters (j) and memory function 
from Eq. 11. As (j) increases, diffusion increases. All curves fit 
closely to quadratic distributions. 



We do this by an appropriate modification of the quan- 
tum walk step: MSC -> MSU{(t))C with the action de- 
fined: 



C 
U 

s 

M 



F,ci, 

|x,ci, 

|x,ci, 
|x,ci, 



3 

^cl]j ((?^) ci , . . . , Cat- 1 , j) , 



,Cn) ^ |x + CAr,Ci, . . . ,CAr), 
,cn) |x,CAr,Ci, . . . ,CAr_i), 



where 



(12) 



(13) 



Here we allow for the rotation angle (j){x) to also be a 
function of position. Note that up to conjugation by a 
local rotation on the A^th coin, this unitary is generated 
by an Ising-like coupling between coins 1, . . . , — 1 and 
the Nth coin. Effectively, the unitary U{(j)) performs a 
conditional rotation on the Nth coin with an angle that 
is proportional to the polarisation 5*^ = cr^ of all 

the other coins in memory. 



VII. QUANTUM WALKS WITH MEMORY ON 
GENERAL GRAPHS 



In our studies so far we have restricted ourselves to 
regular graphs. As discussed, we can generalise to regular 
graphs of constant order by increasing the dimensional- 
ity of the memory states. Next we consider the question 
whether the formalism can be generalised to arbitrary 
graph structures. To address this question we turn to the 
formalism presented in Ref. [27], which introduced an 
approach for modelling multi-walker quantum walks on 
arbitrary graphs. We will restrict ourselves to the case of 



a single walker, in which case the formalism is, 
C: \x,c)^ 



S : \x,j) \j,x). 



(14) 



Here c no longer specifies a direction, but rather a vertex. 
Thus, after the coin operator c can be interpreted as the 
vertex the walker will hop to at the next step, and after 
the step operator as memory of the previous vertex. Ux is 
the neighbourhood of vertex x, i.e. the vertices connected 
to X by an edge. 

Next one might assume we can logically generalise this 
to have memory by letting our state vectors be of the 
form |x,ci, . . . ,CAr), as before. The problem with doing 
this is that the number of possible values of is in 
general larger than the size of the neighbourhood of x, 
\nx\' Thus, no unitary matrix A^^^ can be defined which 
maps the cat vertex memory value to a vertex in the 
neighbourhood of x. Therefore, in this intuitive choice of 
formalism, it is not possible to unitarily define a quantum 
walk with memory. If it is possible at all, a more novel 
approach will have to be employed, or alternately non- 
unitary operators could be employed, which would result 
in mixed states. 

A different approach to considering arbitrary graphs is 
to instead focus on regular graphs, such as a regular 2D 
lattice, but allow the coin matrices A^^^ to be different at 
each vertex. By this approach irregularity can be mod- 
elled, even though the underlying graph is regular. For 
example, to model a vertex which does nothing in a 2D 
lattice, the corresponding coin matrix could be defined 
as the identity matrix A^^^ = I = diag(l, 1, 1, 1) Vx. This 
would have the effect of propagating the walker straight 
through the vertex without entering a superposition. 



VIII. MULTIPLE WALKERS 

Next we generalise our formalism to multiple walk- 
ers. To do this we adopt the walker operator formal- 
ism presented in Ref. [27]. Here we replace position/coin 
state vectors with creation operators (in an optical con- 
text these are photon creation operators), of the form 
iL'(x, ci, . . . , cat)^. Using this approach a system with mul- 
tiple walkers can be easily modelled by adding additional 
walker operators. The transformations on the walker op- 
erators are defined analogously to Eq. 12, 

C : w{x,Ci,...,CNy ^A^^^^^.^(x,Cl,...,CA^-l,j)^ 

j 

S : w{x, ci, . . . , cn)^ w{x + Cat, ci, . . . , CA^)^ 
M : w{x, ci, . . . , Cat)"*" w{x, cat, ci, . . . , cn-i)^ • (15) 

Using this formalism a single walker state would 
be defined as k;(x, ci, . . . , CAr)^|0), where |0) repre- 
sents an empty graph (the vacuum state in an op- 
tical context). With n walkers this generalises to 
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w{x('\c\'\ w{x(-\c^^\. . . , c(^))t|0) (up to 

normalisation), where c'f^ is the coin associated with 
memory element i of walker j. 

We will not present any numerical results for mul- 
tiple walkers here, since the complexity of the system 
grows exponentially with the number of walkers, making 
it computationally intractable for even a modest amount 
of memory. This will be discussed further in Sec. XII. 

Bosons only interfere with one another when they are 
indistinguishable. Thus, in a continuous-time quantum 
walk, where no coins are present, multiple walkers will 
only interfere when they share the same position. In an 
ordinary single-coined discrete-time quantum walk, mul- 
tiple walkers will only interfere when they share the same 
position and coin values. Similarly, multiple walkers with 
memory will only interfere when they share the same po- 
sition and their entire memory history is the same. For 
this reason it is to be expected that there will be much 
less quantum interference taking place between walkers 
with memory, since the probability of different trajecto- 
ries sharing the same memory history drops exponentially 
with the size of the memory. 



IX. MEMORY AS A DECOHERENCE MODEL 

We have seen that full memory reduces our quantum 
walk to a walk with classical statistics. We now develop 
some further intuition as to why this is the case. 

The most general form of discrete-time evolution is 
the quantum process formalism - this formalism captures 
unitary evolution, measurement, as a well as decoherence 
processes. It is well known that any quantum process can 
be expressed as a primary system P interacting with an 
environment system where the environment is traced 
out and only the primary system is observed [ ]. Thus, 
any quantum process can be expressed as, 

p = = trE{Up,E\i')WUl',E)^ (16) 

where Up^e is a unitary acting jointly on the primary 
and environment systems, which in general entangles the 
two systems. When Up^e is an entangling operation, the 
quantum process will in general yield a mixed state. 

In the case of our quantum walk with memory this 
unitary takes the form U p^e = MSG. Here C and M act 
only on the coin subsystem, whereas S is an entangling 
operation that couples the two subsystems. Thus, Up^e is 
entangling. Following evolution we only measure P, trac- 
ing out a model analogous to Eq. 16. This is illustrated 
in Fig. 12. 

To understand how full memory reduces a quantum 
walk to have classical statistics, let us consider the first 
step in the evolution of a quantum walk. Since we are just 
considering the first step, it suffices to introduce just one 
coin, in which case the evolution over one step is given 



position 



CM 




measurement 

e((M5C)*|^in>(^in|(Ct5tMt)*) 



|ci, . . . ,Cn)^ 

memory 



FIG. 12: Visualising memory as an environment system. The 
coin and memory update operators act only on the memory 
subsystem, while the step operator couples the memory en- 
vironment with the position. At the end of the evolution we 
trace out the environment, only observing the position sub- 
system, yielding, in general, a mixed state. 



by, 



I^Ain) = |0,+1), 

C'li'in) = (|0,+l)-i|0,-l))/V2, 
5C|Vi„) = (|l,+l)-i|-l,-l))/\/2. 



(17) 



After SC we have a maximally entangled state across 
two positions. Tracing out the coin system we obtain the 
maximally mixed state across two positions, 

tr^(5C|^in){^i„|Ct5t) = (|1){1| + I - l){-l|)/2, (18) 

what one would obtain for a single step of a classical ran- 
dom walk. Obviously if we repeat this many times, each 
time employing a fresh coin, our evolution will proceed 
classically, since introducing a fresh coin, maximally en- 
tangling it with the position and then tracing out the 
coin is simply a balanced classical coin-flip. 

Thus, full memory, which is equivalent to introducing 
a new coin at each step, yields a classical walk, provided 
the coins are all traced out prior to measurement of the 
position. On the other hand, with partial memory we 
expect to see purely classical statistics for t < after 
which quantum behaviour emerges since there are quan- 
tum correlations between past and present coin values. 
With no memory (i.e. a single coin) we are re-coupling to 
the same environmental degree of freedom at each step, 
and thus we do not reduce to a classical walk since each 
coin-flip is not independently random. 

From a non-rigorous, intuitive point of view, this re- 
sult is perhaps not surprising. As we increase the size 
of the environment we are effectively making the joint 
system more macroscopic, and increased decoherence is 
intuitively not surprising. 

Based on the arguments presented, we expect that de- 
coherence is a generic feature of quantum systems with 
memory where the memory is not measured. 
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Time 

FIG. 13: (Color online) Comparison of various quantum walk 
behaviours. Quantum walks on a line with: five recycled coins, 
memory function from Eq. 9, no randomness (long dashed); 
one coin and no randomness, i.e. the standard quantum walk 
(solid); five recycled coins, trivial memory function, and no 
randomness (short dashed); five recycled coins, memory func- 
tion Eq. 9, spatially random coin flip (dotted); five recy- 
cled coins, trivial memory function, spatially random coin 
flip (dot-dashed); one coin, spatially random coin flip (dot- 
dot-dashed). The plots with randomness were generated by 
averaging over 30 realisations of random spatially dependent 
angles for the coin flip, where the mean is plotted in the mid- 
dle and one standard deviation plotted above and below. 

X. WALKS WITH SPATIAL RANDOMNESS 

It is known that spatial randomness in various forms 
can give rise to localization of quantum walkers. In par- 
ticular, using a spatially dependent coin with rotation 
angles that are chosen randomly and quenched, leads to 
a localized distribution. Localization is a phenomenon 
that arises due to coherence, specifically destructive in- 
terference between amplitudes for propagation away from 
the origin. Joye and Merkli [ ] and Ahlbrecht et al. [ 5] 
show that dynamical localisation occurs for a spatially in- 
homogeneous coin, when the coin parameters are chosen 
randomly from continuous or discrete sets. If temporal 
randomness is also present, e.g. when the coin parame- 
ters change in time, then localisation does not occur and 
instead the quantum walk looks classically diffusive [ ] . 
In contrast, the walker spreads ballistically or diffusively, 
but no localisation has been observed. 

We examine how the presence of a memory affects such 
behaviour. To do so we use the memory function in Eq. 
9 but with an angle that varies randomly in space: 

^=i+^(-)i^7f^ (19) 

where 4>{x) is chosen randomly in the range 
— 1 < (j){x) < 1, independently at each position. 

Results for various types of walks are shown in Fig. 13. 
The three walk scenarios without spatial randomness all 
behave quantum mechanically, i.e. they have quadratic 



dispersion. The walker with 5 coins and a memory func- 
tion spreads notably faster than that without, and in- 
deed after a crossover time of t = 15 spreads faster than 
the standard quantum walk with one coin. We cannot 
confirm that this is universal behaviour, however for var- 
ious initial conditions the walk with memory function 
does spread quadratically, and it was already proven that 
the walk with finite number of spins (less than t) and 
trivial memory function also spread quadratically [38]. 
The results with spatial randomness are more striking. 
The one coin walk with randomness shows localisation 
as predicted. However the 5 coin walk with trivial mem- 
ory function does not localise and indeed has a variance 
scaling linearly with t. We interpret this as due to the 
fact that the multiple coins provide a sufficiently large 
environment such that the coherence necessary to localise 
is lost. This is in contrast to the case without random- 
ness where enough coherence is preserved to maintain 
quadratic dispersion as discussed above. In effect this is 
indicating that the coherence necessary to observe locali- 
sation is more sensitive to decoherence than is the coher- 
ence to obtain quantum mechanical speedup. Even more 
surprising is that the quantum walk with memory func- 
tion spreads much faster in the presence of randomness 
than either of the other two cases. For this size simula- 
tion we are unable to confirm the asymptotic behaviour 
of this case though the linear coefficient of the variance 
does dominate for t < 80. We conjecture the sub-ballistic 
behaviour of the multi-coined walks in the presence of 
randomness as a consequence of the coins introducing an 
effective temporal randomness to the dynamics by rein- 
troducing past spatial randomness at each new step of 
the walk. The memory function we have used tends to 
average out the effect of spatial randomness from earlier 
steps which could explain why it is spreading faster than 
the case of the walk with trivial memory function. 



XL EXPERIMENTAL CONSTRUCTION 

Linear optics quantum computing (LOQC) [47-49] has 
become one of the more promising candidates for im- 
plementing large-scale quantum information processing 
devices. We now briefly discuss the prospects for linear 
optics implementation of quantum walks with memory. 

When using single photons to represent walkers, this 
formalism can always be experimentally constructed us- 
ing just passive linear optics networks (i.e. beamsplitters 
and phase-shifters) in a manner similar to that described 
in Ref. [ ]. As a simple example, the explicit linear op- 
tics construction of an TV = 2 quantum walk with two 
positions is depicted in Fig. 14. Each 'wire' represents an 
optical mode, and wires are bundled together to repre- 
sent the different position/coin combinations. In Fig. 14 
the coin operator acts only on the C2 space, while S and 
M are permutation operators as per Eq. 15. 

In general, with an A/'-element memory, d position 
states and |c|-dimensional coins, there will be d - \c\^ 
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wires. Thus the number of modes required for experimen- 
tal implementation grows exponentially with the length 
of the memory, making large- memory demonstrations un- 
viable. 




FIG. 14: Explicit linear optics construction of an N — 2 quan- 
tum walk with two position states. The 'wires' represent op- 
tical modes and are bundled together to represent different 
configurations of position and coin values. The coin operators 
act on C2 only, while S and M are permutation operators. 

Experimental construction of wise of man walks would 
be particularly challenging, as it would require imple- 
menting conditional coin operations, coherently con- 
trolled by some arbitrary function of the memory history. 
This is largely beyond the scope of present-day linear op- 
tics implementations. 

Nonetheless, linear optics demonstrations of walks with 
small memories and without history-dependence are vi- 
able in the near future. An experiment like that shown 
in Fig. 14 could be readily constructed with present-day 
technology. 



XII. CLASSICAL SIMULATION & 
COMPUTATIONAL COMPLEXITY 

As discussed earlier, the number of optical modes re- 
quired to implement a quantum walk with memory grows 
as d- |c|^. This makes experimental construction chal- 
lenging, since the number of optical modes grows expo- 
nentially with the size of the memory. With n walkers, 
the number of unique basis states grows as 0(((i • |c|^)^). 
Thus, the complexity of brute- force classical simulation 
of this quantum system grows exponentially against both 
the number of walkers and the length of the memory. 
Therefore classical simulation of such systems will be lim- 
ited to walks with a modest number of walkers and small 
memories, and for this reason we do not present numer- 
ical results for mult i- walker scenarios. Results for multi- 
ple walkers without memory were previously presented in 



Ref. ["T. 

From an algorithmic perspective, we are interested in 
systems with exponential complexity, since otherwise the 
system will be classically simulatable and exponential 
quantum speedup will obviously not be possible. It was 
noted in Ref. [ ] that if the size of the graph is 'effi- 
cient' (i.e. d- \c\^ = 0(poly(p)), where p is the size of 
the problem), then the only way to achieve exponential 
complexity is with the addition of multiple walkers, as 
noted above. 

Because quantum walks with memory can be imple- 
mented with just linear optics, then, with multiple walk- 
ers, such systems are a subset of the Boson-sampling 
problem described by Aaronson & Arkhipov [ ], who 
presented strong evidence that such systems are classi- 
cally hard to simulate. It was shown in Ref. [50] that 
an ordinary single-coined quantum walk with multiple 
walkers is already universal for Boson-sampling (com- 
plexity class denoted BosonSampP). Thus, at best, 
adding memory to the walk will also be universal for 
Boson-sampling - adding memory cannot extend the 
complexity class of the system beyond BosonSampP. 
It is known that BosonSampP C SampBQP (the sam- 
pling version of BQP), and it is strongly believed that 
BosonSampP C SampBQP [ ]. Therefore we expect 
multi-walker quantum walks with memory will present 
no algorithmic advantage over multi- walker single-coined 
quantum walks, although they are nonetheless likely 
to be classically hard to simulate, and quantum walks 
with memory (on efficiently sized graphs) are likely not 
universal for quantum computation. However, on expo- 
nentially large graphs it is known that both discrete- 
and continuous-time quantum walks without memory are 
BQP [17, 18], and we don't rule out that the same ap- 
plies to quantum walks with memory. 



XIII. CONCLUSION 

We have presented a formalism for quantum walks 
where the walker has arbitrary amounts of memory of 
its previous history. We examined the effects of history- 
dependent coins and introduced the memory function, 
which manipulates the diffusive properties of the walk 
as a function of the walker's history. With appropriate 
choice of memory function, diffusion can be reduced, en- 
hanced, or be made symmetric or asymmetric. We also 
studied the effects of measurement of the memory reg- 
isters on the evolution of the walk. Measurement of the 
oldest memory elements has the effect of 'resetting' the 
quantum walk to begin at a new position that is a func- 
tion of the measurement outcomes. In the case of a two- 
dimensional walk on a regular rectangular lattice, mem- 
ory eliminates spatial entanglement between the two spa- 
tial degrees of freedom, yielding a separable spatial dis- 
tribution, even if entangling coins are employed. 

A key observation from our numerics arises in the con- 
text of spatial randomness. We found that in settings 
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where single coined walks localize, multi-coin walks do 
not. Moreover, walks with non-trivial memory function 
appear to spread faster than those which only recycle 
coins. We conjecture that this occurs due to the the 
memory function averaging out the mechanism which in- 
duces an effective temporal randomness to the walk due 
to recycling coins. It has been shown that for one dimen- 
sional walks with both spatial and temporal randomness, 
the former dominates to lead to classical diffusion [5c]. 
Whether walks with memory can counteract this effect is 
unknown. 

Finally, we explicitly demonstrated how to physically 
construct an optical quantum walk with memory. How- 
ever the required physical resources grow exponentially 
with the length of the memory. Thus, experimental con- 
struction of a quantum walk with large memories is lim- 
ited. We presented a formalism for mult i- walker quan- 
tum walks with memory and discussed the prospects for 
classical simulation. Unfortunately the required classical 
resources grow exponentially against both the number of 



walkers and the size of the memory. Therefore classical 
simulation is challenging beyond a few walkers with even 
modest memory. We argued that a multi- walker quantum 
walk with memory has no algorithmic advantage over its 
single-coined counterpart. Algorithmic complexity aside, 
quantum walks with memory may be interesting to pur- 
sue nonetheless owing to their unique physical properties, 
such as their rich diffusive characteristics, which may be 
interesting for quantum simulation applications. 
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